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ABSTRACT 

We carry out a set of self-consistent TV-body calculations to compare the decay rates of 
satellite dwarf galaxies orbiting a disc galaxy embedded in a dark matter halo (DMH). 
We consider both spherical and oblate axisymmetric DMHs of aspect ratio qt = 0.6. 
The satellites arc given different initial orbital inclinations, orbital periods and mass. 
The live flattened DMHs with embedded discs and bulges are set-up using a new fast 
algorithm, MaGalie (Boily, Kroupa & Penarrubia 2001). 

We find that the range of survival times of satellites within a flattened DMH 
becomes ~ 100% larger than the same satellites within a spherical DMH. In the oblate 
DMH, satellites on polar orbits have the longest survival time, whereas satellites on 
coplanar prograde orbits are destroyed most rapidly. The orbital plane of a satellite 
tilts as a result of anisotropic dynamical friction, causing the satellite's orbit to align 
with the plane of symmetry of the DMH. Polar orbits are not subjected to alignment. 
Therefore the decay of a satellites in an axisymmetric DMH may provide a natural 
explanation for the observed lack of satellites within — 30° of their host galaxy's disc 
(Holmberg 1969; Zaritsky & Gonzalez 1999). 

The computations furthermore indicate that the evolution of the orbital eccen- 
tricity e is highly dependent of its initial value e(t — 0) and the DMH's shape. We 
also discuss some implications of flattened DMHs for satellite debris streams. 

Key words: stellar dynamics - methods: numerical - galaxies: evolution - galaxies: 
kinematics and dynamics - galaxies: spiral - dwarf satellites 



1 INTRODUCTION 

Non-spherical mass distributions around galaxies and galaxy 
clusters are needed to reconcile the dynamics with lensing 
statistics (Mailer et al. 2000; Gonzalez et al. 1999; Mailer 
et al. 1997; Keeton & Kochanek 1998) or galactic disc 
warps (Binney 1992). In CDM cosmogony, aspherical bound 
dark matter haloes (= DMHs) form as a result of gravi- 
tational clustering. Dubinsky (1994) finds in his computer 
simulations a Gaussian distribution of DMH aspect ratios, 
qh = c/a > 0, where c and a are the minor and major axes 
of an oblate spheroid, of mean < qh >= 1/2 and dispersion 
equal to 0.15. In a CDM framework, therefore, DMHs may 
achieve an aspect ratio as high as qh = 0.65, covering a range 
of values more than adequate to account for lensing data: 
for instance, Mailer et al. (2000) find that for the galaxy 
B1600 + 434, 0.5 < q h < 0.75. The inferred morphology of 
the DMH, however, depends on the details of its radial mass 
profile and embedded baryonic galaxy components, such as 
the disc and bulge of a spiral galaxy, and is, on the whole, 
a loosely constrained quantity. For instance, Oiling & Mer- 
rifield (2000) studying the axis-ratio of the Milky Way by 



two independent methods (first, by the measurements of its 
rotational curve and the amount of dark matter in the solar 
neighborhood and, second, by the variation in thickness of 
HI emission) find a consistent value of qh ~ 0.8, although it 
depends strongly on the measurement of the Milky Way's 
parameters. 

In general, different techniques yield a wide range of 
values for qh- Models of decaying neutrinos (Sciama 1990) 
or cold molecular gas (Pfenniger et al. 1994) suggest a mini- 
mum value qh ~ 0.2, whereas models of the Milky Way halo 
suggest qn — 0.9 to account for the narrow debris stream 
of the Sagittarius dwarf (Ibata et al. 2001). The measure of 
axis-ratios throughout a DMH is riddled with uncertainties 
due to the unknown profile of the halo. Furthermore, lensing 
or stellar kinematics yield diagnostics that are constrained 
in the inner region mainly, where density gradients are the 
largest. The Sagittarius dwarf data suggest that the inner 
region of the MW DMH would be spherical, yet questions 
arise as to whether the same hold true on scales approaching 
100 kpc or more. 

The dynamics of galactic satellites may help constrain 
DMH profiles on large scales by direct observations of their 
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distribution around the host galaxy. Holmberg (1969) and 
Zaritsky & Gonzalez (1999) point out that satellites around 
disc galaxies are found more often aligned with the poles of 
the host galaxy, the so-called 'Holmberg effect'. One possible 
reason for this effect is the enhanced satellite-disc coupling 
for co-planar satellite orbits. Quinn & Goodmann (1986), 
however, find in their iV-body study that discs alone cannot 
account for the original statistical distribution of Holmberg's 
data. 

A remedy may be sought in the form of an extended 
non-spherical DMH. An anisotropic velocity (and mass) dis- 
tribution will cause a satellite's orbit to align with the axes 
of the velocity ellipsoid of the host galaxy (Binney 1977). 
A strategy for exploring the parameter space of orbits of a 
population of satellites would be to integrate the orbit of 
satellites in a fixed non-spherical potential using e.g. Bin- 
ney's treatment for the dynamics of a point source. In the 
last years, there have been several studies of how dynamical 
friction influences the orbit and structure of satellites that go 
in this direction (e.g. Colpi et al 1999 using the theory of Lin- 
ear Response), as well as N-body simulations (e.g. van den 
Bosch et.al. 1999). However there is at present no body of 
work for this problem against which to compare the analytic 
treatment, as done for spherically symmetric systems (see for 
instance Taylor & Babul 2000). Furthermore, mass loss by 
the satellite is difficult to account for analytically (see Helmi 
& White 1999; Johnston et al. 1999). In addition, fixed po- 
tentials prevent important feed-backs from the dynamical 
friction process. For example, Weinberg (2000) argues that 
the wake from a heavy satellite may induce bending modes 
in the disc. Effects such as wake, disc bending, motion of 
the primary galaxy's centre of mass, would in turn influence 
the satellite's orbit. 

In this paper, we study how axisymmetric (flattened) 
DMHs affect the orbital decay and survival of satellites, 
paying particular attention to the orbital inclination of the 
satellite with respect to the disc and DMH plane of symme- 
try. We are motivated by the Holmberg effect, and by the 
fact that no study of satellite decay exists to date which 
takes into account both the velocity anisotropy and the flat- 
tened density structure of the DMH. Given the likely im- 
portance of feedback, and the as yet untested analytical de- 
scription of dynamical friction in anisotropic systems, we 
resort to fully self-consistent calculations with live multi- 
component galaxies and satellites. This has become possible 
only very recently through the availability of a new algo- 
rithm, MaGalie, that allows the construction of large-N 
flattened multi-component galaxy models (Boily, Kroupa & 
Penarrubia 2001, hereinafter BKP). 

Section 2 introduces the models. In Section 3 we study 
how flattened DMHs affect satellite decay, especially in com- 
parison with spherical DMHs without velocity anisotropy. 
We also touch on orbital precession and its implications for 
the spread of tidal debris. The paper concludes with Sec- 
tion 4. 



2 THE MODELS 

A subset of our spherical models are similar to the models 
of Velazquez & White (1999, hereinafter VW) to facilitate 
an inter-comparison of different numerical treatments. 



2.1 The primary galaxy model 

In order to minimize computational time when constructing 
flattened DMHs with embedded bulges and discs, we apply 
a new highly-efficient technique using multi-pole potential 
expansions to tailor the local velocity ellipsoid to the re- 
quired morphology (BKP). The algorithm to add together 
individual components in a single galaxy is adapted from 
Hernquist's method (Hernquist 1993). The new code, Ma- 
Galie, scales linearly with particle number and hence we 
can construct flattened DMHs consisting of <; 10 6 particles 
or more, in a short computational time. 

For the density distributions of the disc we take 

p d (R,z) = — — exp(-R/R d )sech 2 (z/z ), (1) 

Md being the disc mass, zq the vertical thickness, and Rd 
the exponential scale length in the radial direction. The 
mass profile decays exponentially with R and is composed of 
isothermal sheets along the vertical direction. Velocities are 
assumed to have a GauBian distribution. The square of the 
radial velocity dispersion is taken to be proportional to the 
surface density (see Lewis & Freeman 1989), v'^ oc E(i?) = 
E(0)exp(— R/Rd), where the constant of proportionality is 
determined by fixing Toomre's Q-parameter at the Solar ra- 
dius. Following VW, we select Qq = Q(R@) = 1.5. The ver- 
tical component of the velocity ellipsoid is v 2 = nCE(R)zo 
in agreement with an isothermal sheet (Spitzer 1942). The 
azimuthal component is obtained from the epicyclic approx- 
imation cr? = v^k 2 / (4Q. 2 ) (e.g. Binney & Tremaine 1987, 
hereinafter BT). 

For the bulge we adopt the spherical Hernquist profile 
(Hernquist 1990), 

_M b a 

Pb ~ 2tt r(r + a)3' [Z> 

where Mj is the bulge mass and a the spherical scale length. 
This analytical profile fits the de Vaucouleurs law (de Vau- 
couleurs 1948). The velocity field is constructed from the 
Jeans equations by assuming isotropic Gaussian velocity dis- 
tributions at each radial distance (Hernquist 1993). 

We use a non-singular isothermal profile for the DMH, 

= M h a exp(-r 2 /r 2 ut ) 
27r 3 / 2 r cu t r 2 + 7 2 
M h being the DMH mass, 
core radius, and 



(3) 

a t the cut-off radius and 7 the 



-2)f3 2 + Q{f3 3 ) 



(4) 



a = {l- v^73exp(/3 2 )[l 
1 + + (tt 

where f3 = 7/r cut 1/24 in our calculations. For j3 = 1/24 
we find a ~ 1.076 — > 1 already and hence thereafter we set 
q = 1 in our analysis. To construct the flattened (oblate) 
DMHs, a non- homologous transformation is applied to (Bl) 
to achieve the desired axis ratio c/a while preserving the 
central density: 1) first we flatten the DMH down the z- 
axis only, until the desired aspect ratio qu — c/a is reached; 
2) then all axes are stretched by a factor A, A = XR, such 
that M/R 3 = M/A 2 C = M/[A 3 q h ] = M/[R 3 \ 3 q h ]. Solving 

— 1/3 

for X — q h m 1.18 when qu = 0.6. The orbital period 
oc 1/y/Gp, and hence the dynamical time, is unchanged in 
the new equilibrium. Particle velocities are obtained by ad- 
justing the initial isotropic distribution (as for the bulge) to 
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Figure 1. Total contribution from the three Gl galaxy compo- 
nents (disc, bulge and halo) to the circular velocity (solid line). 
We also plot the circular velocity for each galaxy component. On 
very small scales (r < 1 kpc) the bulge accounts for the bulk 
of V c . Further out, the dynamics is dominated by the halo. The 
solar radius is at Rq = 8.5 kpc. 

the oblate iso-potential surfaces leading to a stable axisym- 
metric DMH with embedded bulge and disc (BKP). 

We use four different isothermal DMH models: a spher- 
ical (Gl) and a flattened DMH (G2) with axis-ratio qu = 
c/a — 0.6, which lies within the distribution of Battenings 
given by CDM models. The third (G3) and fourth (G4) mod- 
els have the same properties as Gl and G2, respectively, but 
with enlarged cut-off radii. We can define two typical dis- 
tances, the core radius in the symmetry plane (7 a ) and in 
the vertical direction (-y c ). Since concentric iso-density con- 
tours have the same axis-ratio throughout the DMH, both 
core radii are related by qn = Ic/la- 

Our system of units is such that Md = Rd = 1 and 
G = 1. According to Bahcall, Smith & Soneira (1982), 
M d = 5.6 x 1O 1O M and R d = 3.5 kpc for the Milky Way 
which we adopt as a typical primary galaxy model, so that 
time and velocity units are, respectively, 1.3 x 10 7 yr and 
262 kms -1 . The half-mass radius of the disc is located at 
R0.5 ~ 1.7 Rd = 5.95 kpc, with a rotation period of 13 time 
units. Table 1 summarizes the parameters and Fig. ^ plots 
rotational curves for two models. 

As VW point out, there are some caveats to keep in 
mind concerning the above models: (i) The DMHs are pos- 
sibly too small in mass and extension. Zaritsky & White 
(1994) show, by studying satellite orbits in the Local Group 
and external galaxies, that DMH limits may extend beyond 
200 kpc with masses over 2 x 1O 12 M0. However, as VW 
comment, the velocity curves of our DMHs Gl and G2 are 
consistent with the largest velocities observed for stars in 
the solar neighbourhood (Carney & Lathman 1987), and 
they are possibly massive enough to give realistic velocities 
of satellites on eccentric orbits, (ii) The DMHs may be too 
concentrated. Persic, Salucci & Stel (1996) argue for a DMH 
core radius of 7 = (1 — ► 2) X ir. op t, R op t = 3.2 Rd , where Rd 
is the disc scale-length. However the DMH parameters were 
selected to avoid bar formation in the disc. We observed 
that a less concentrated DMH or bulge allows a stable disc 
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Symbol 


Value(ph.u) 


Value (m.u) 


Disc 




100000 






M d 


5.60 x 10 10 M Q 


1.00 




Rd 


3.50 kpc 


1.00 




zo 


1.40 kpc 


0.40 




y© 


l.OU 


l.OU 




Fir 




2.43 


Bulge 


N b 


33328 






M b 


1.87 x 10 10 M Q 


1/3 




a 


0.53 kpc 


0.15 


DMH (Gl) 


N h 


1400000 




(spherical) 










M h 


7.84 x 10 n M Q 


14.00 




1 


3.50 kpc 


1.00 




1h 


1.00 


1.00 




rcut 


84.00 kpc 


24.00 


DMH (d2) 


Nu 


1400000 




(oblate) 










M h 


7.84 x 10 n M Q 


14.00 




Qh 


0.60 


0.60 




la 


3.80 kpc 


1.10 




7c 


2.28 kpc 


0.65 




Tcut 


84.00 kpc 


24.00 


DMH (G3) 


N h 


1400000 




(spherical) 










Mh 


7.84 x 10 11 Mq 


14.00 




7 


3.50 kpc 


1.00 




Qh 


1.00 


1.00 




r cu t 


133.00 kpc 


38.00 


DMH (G4) 


N h 


1400000 




(oblate) 










M h 


7.84 x 10 n M Q 


14.00 




Qh 


0.60 


0.60 




la 


3.80 kpc 


1.10 




7c 


2.28 kpc 


0.65 




fcut 


133.00 kpc 


38.00 



Table 1. Primary galaxy models. Oblate models have an aspect 
ratio = 0.6. The units are such that Ph.u. means 'physical 
units', and m.u. 'model units'. 

to form a bar after few satellite passages. With our 7, the 
presence of a bar is avoided at least until the destruction of 
the satellite. 



2.2 Satellite models 

We use self-consistent King models (King 1966) to represent 
our dwarf galaxies. These models fit early-type dwarf galax- 
ies (Binggeli et al. 1984), where r c and r t are the core and 
tidal radii, respectively. For a comparison with the work of 
VW we adopt c = 0.8. 

To construct the models we choose the satellite mass 
M s , r c , rt and thus c. The tidal radius is determined by 
computing the density contrast, Ps(rt)/~fhj(r a ) ~ 3, at the 
apo-centric distance (r a = 55 kpc) at t = 0, Ps( r ) be- 
ing the averaged density of the galaxy (same procedure as 
VW). This guarantees that all satellite particles are bound 
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Symbol Value(ph.u) Value (m.u) 




iV s 


40000 




M a 


5.60 x 10 9 M Q 


0.10 


*(0)/o-2 


5.00 


5.00 


r c 


1.00 kpc 


0.29 


Tt 


6.31 kpc 


1.80 


c 


0.80 


0.80 


< r > 


1.64 kpc 


0.47 


erg 


52.00kms~ 1 


0.20 


iV s 


40000 






1.12 x 1O 1O M 


0.20 


*(o)/«g 


5.00 


5.00 




1.00 kpc 


0.29 


n 


6.31 kpc 


1.80 


c 


0.80 


0.80 


< r > 


1.64 kpc 


0.47 


O"0 


74.00kms~ 1 


0.28 



Table 2. Satellite models. *(0) = *(r t ) - *(0), *(0) being the 
central potential and < I > (rt) the potential at the tidal radius (fol- 
lowing BT notation); tro is the velocity dispersion at the centre, 
and < r > the average radius of the satellite. 



2 4 8 

r (kpc) 

Figure 2. Rotational curve of the satellite models SI and S2 (see 
Table 2 for the characteristics of each one). 



at t = 0. Tables for the numerical rendition of the corre- 
sponding King profiles can found in BT or in the original 
paper of King (1966). Table 2 summarizes the parameters, 
while Fig. [| plots rotational curves. Note that we use the 
same M s , r c and "rt" despite placing t he s atellites at differ- 
ent apogalactica r a > 55 kpc (Section 2.3), which increases 
the true tidal radius of the satellite, though the stability 
condition at t = is still well- accomplished. We do this 
rather than using different r c or r t in order to study the 
same satellites on different orbits. 

Our satellites are much more massive than the Milky 
Way dSph satellites which have M s < 10 8 M , but our 
adopted values are typical for the satellites that enter dis- 
tant samples such as used by Holmberg (1969) and Zaritsky 
& Gonzalez (1999). 



2.3 Numerical method and orbital parameters 

We use Superbox (Fellhauer et al. 2000) to evolve the 
galaxy-satellite system. Superbox is a highly efficient par- 
ticle mesh-code based on a leap-frog scheme, and has been 
already implemented in an extensive study of satellite dis- 
ruption by Kroupa (1997) and Klessen & Kroupa (1998). 

Our integration time step is 0.39 Myr which is about 
l/25th the dynamical time of satellite S2. We have three 
resolution zones, each with 64 3 grid-cells: (i) The inner grid 
covers out to 3 radial disc scale-lengths, which contains ss 
90 % of the disc mass, providing a resolution of 350 pc per 
grid-cell, (ii) The middle grid covers the whole galaxy, with 
an extension of 24 disc scale-lengths (84 kpc) for the models 
Gl and G2, giving a resolution of 2.8 kpc per grid-cell. The 
satellite always orbits within this grid except when it reaches 
the disc, avoiding cross-border effects (see Fellhauer et al. 
2000). For the models G3 and G4, the middle grid extends 
to 141 kpc and has a resolution of 4.7 kpc per grid-cell. The 
orbits of the satellites are located within this zone, (iii) The 



outermost grid extends to 348 kpc and contains the local 
universe, at a resolution of 11.6 Kpc. 

As for the satellite grid-structure, the resolutions are 
816 pc per grid-cell for the inner grid that extends to 
24.48 kpc, 1.2 kpc per grid-cell for the middle grid which 
extends to 36 kpc, and 11.6 kpc per grid-cell for the outer- 
most grid that covers the local universe. Only the inner and 
middle grids move along with the satellites, remaining po- 
sitioned on their centre-of-density locations. The outer grid 
is identical for primary galaxy and satellite. 

Klessen & Kroupa (1998) compared calculations per- 
formed with SUPERBOX with direct-integration N-body 
calculations and found good agreement. Specifically, they 
verified that varying the grid resolution by factors of a few 
did not lead to unstable satellite models. The stability of 
the satellite models does not depend strongly on the val- 
ues adopted here. Furthermore, based on the comparison 
with the direct-integration method, the heating introduced 
by two-body effects prove entirely negligible for the model 
satellites we consider. The selection of grid parameters en- 
sures the conservation of energy and angular momentum for 
satellites in isolation over times as long as our calculations 
to a high degree. Conservation of total energy and angular 
momentum is better than 1% for all the models. 

The disc is poorly resolved in the z-direction and we 
do not study its evolution in any detail. We verified that 
the disc parameters do not evolve for galaxies in isolation 
(no satellites). Since Superbox is a mesh code, a poor z- 
resolution for the disc is expected due to the limited num- 
ber of grids. This provokes the disc modeled here to be un- 
realistically thick, however it does provide a quadrupolar 
(non-spherical) potential of the appropriate magnitude. A 
mesh code has the advantage that it does not introduce self- 
heating since it does not calculate two-body interactions, 
which would have been significant in the disc given the fi- 
nite number of particles used (see the discussion in VW). 
The effects on the satellite dynamics due to two-body in- 
teractions are drastically reduced by the low mass of the 
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halo particles (see Steinmetz & White 1997). Furthermore 
the disc heating by halo particles is minimized since each 
component particle masses are in a one-to-one ratio. 

The decay of satellites with various masses through dy- 
namical friction in an extended spherical DMH is studied by 
Fellhauer et al. (2000), who found good agreement between 
SUPERBOX calculations and Chandrasekhar's (1960) for- 
mula for dynamical friction when the Coulomb logarithm A 
is set to In A = 1.6 (cf. eq [| below). This agrees with the 
findings by VW for similar calculations. 

We carry out a set of calculations varying the parame- 
ters of the satellite and the primary galaxy that influence the 
satellite-primary galaxy interaction. These parameters are: 
(i) the initial orbital inclination (#;), defined as the angle 
between the initial angular momentum vector of the satel- 
lite and the initial angular momentum of the disc, (ii) the 
satellite's mass, (iii) the satellite's apo-galactic distance, (iv) 
its orbital eccentricity, and (v) the DMHs ellipticity, 1 — qn- 

Before injecting the satellite into the primary galaxy 
we allow the galaxy and satellite to settle into a stationary 
state by integrating the isolated systems for a few dynamical 
times with Superbox (as in Kroupa 1997). Examples of the 
stationarity of multi-component galaxies are given in BKP. 
The satellite is then placed at apo-galacticon with a velocity 
as described next. 

The orbit of the satellites are rosettes. VW define the 
'circularity' of the orbit as ej = J/ Jc{E), J being the satel- 
lite's angular momentum and Jc the corresponding angular 
momentum for a circular orbit with the same energy E as 
the satellite's orbit. In practice, we take the circular velocity 
from the rotational curve plotted in Fig. [j] at the satellite's 
initial distance, and multiply it by ej. This procedure gives 
an eccentric orbit with the same energy. The parameters of 
the numerical experiments are listed in Table ^. 



3 SATELLITE DECAY 

We discuss our results in general terms below before go- 
ing into detailed consideration of the mass loss and survival 
of satellites (Section 3.2), and the orbital evolution of the 
inclination angle, eccentricity and precession, respectively 
(Sections 3.3 to 3.5). Section 3.6 takes a brief look at the 
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G1S145 
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implications for tidal streams of dissolving satellites. 



Table 3. The numerical experiments. The peri- and apo-galactica 
are r p and r a , respectively, and e = 1 — r p /r a is the orbital 
ellipticity (BT, p.21). 



3.1 Introductory comments 

We denote by 'G1S145' the compound primary galaxy made, 
in this case, of a spherical DMH plus embedded disc and 
bulge, Gl, and satellite SI, in an orbital plane initially set 
at an inclination angle 9 = 45° with respect to the plane of 
symmetry of the system. In what follows we take this model 
as reference, but all models followed a similar evolution. 

There are two main physical mechanisms that regulate 
the satellite's orbital decay: (i) dynamical friction from the 
disc, bulge and DMH, and (ii) tidal interactions, causing in- 
ternal heating and mass loss. The evolution of the satellite's 
orbital radius and mass profile highlight the basic charac- 
teristics of these two processes. Dynamical friction causes 
a steady decrease of the satellite's apo- and peri-centres in 
time as shown on Fig. 3 (dotted line). (Lengths are given in 
model units on the figure but the time is in Gyr.) From t = 



and until t < 2Gyr, both quantities, apo- and peri-centres, 
decrease monotonically. When t > 2 Gyr, the orbital radius 
r ~ 5 or smaller, and the orbital decay is not monotonic 
anymore. The proximity to the disc means that non-radial 
forces affect strongly the remaining evolution, along with 
the structure of the satellite. 

To measure changes in the structure of the satellite, we 
plotted the ten-percentile Lagrange radii centred on the den- 
sity maximum of the satellite (Fig. 3, solid lines). At t ~ 3 
Gyr, the galactic tidal field has inflated the satellite to the 
extent that half of its initial mass is spread throughout the 
volume circumscribed by its orbit. We note that the inner 
10 per cent Lagrange radius is largely unaffected until the 
very late stages of integration. Our strategy for determining 
the orbital parameters of the satellite therefore consisted in 
locating the position of the density maximum of the inner- 
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Figure 3. Evolution of the satellite's Lagrange radii (solid curves, 
defined as the radius at which the spherically enclosed mass 
amounts to 10%, 20%...., 90%) for the model G1S145. The dotted 
line represents the distance of the satellite's centre of density to 
centre of the primary galaxy. Distances are in model units. The 
overall evolution is similar in all other models (Table ^) . 



most Lagrange radius, which then defines a reference coor- 
dinate. 



3.2 Mass loss and disruption times 

To calculate the mass remaining bound to the satellite, 
M s (t), we compute the potential energy <I>i < of each satel- 
lite particle presumed bound to the satellite, and its kinetic 
energy (Ti) in the satellite frame. Following VW, particles 
with Ei = Ti+m s (&i+& cx t) > are labeled unbound, where 
m 3 is the mass of one satellite particle. Particles with Ei > 
are removed and the procedure repeated until only negative 
energy particles are left. $ cxt = GM g (r < r s )/r s > is 
the external potential from the primary galaxy at the satel- 
lite's centre-of-density (r s ). All the particles of the satellite 
are thus assumed to feel the same external potential, which 
is a useful and sufficiently accurate approximation, taking 
into account that most of the bound particles are located 
very close to this point. For example, in Fig. ^| most of the 
satellite's mass lies at a distance less than 4 kpc from the 
position of the centre-of-density until the satellite's disrup- 
tion. This approximation fails whenever the satellite's size 
is comparable to its distance to the galaxy centre. 

Satellites lose mass due to the galaxy's tidal forces. The 
mass loss happens mostly at perigalacticon, since the gradi- 
ent of the galaxy's gravitational force reaches a maximum 
at that point (see Fig. This is seen indirectly in the os- 
cillations of Lagrange radii, always in phase with the orbit 
of the satellite: the satellite fills its Roche lobe and conse- 
quently responds strongly to the changing tidal field. Thus a 
decrease of the apo-galacticon distance implies an enhanced 
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Figure 4. a: Evolution of the satellite mass for M s = O.lMj and 
eccentricity e ~ 0.7. 
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Figure 4 - continued b: As Fig. Ma for satellites with M s = 
0.1 and initial eccentricity e ~ 0.45 and e = 0. (Note that the 
time-axis has changed scale.) 



mass loss. The evolution of satellites exposed to strongly 
varying tidal fields is discussed at length by Piatek & Pryor 
(1995) for one perigalactic passage, whereas long-term satel- 
lite harassment is addressed by Kroupa (1997) and Klessen 
& Kroupa (1998). Consequently, we will not study the in- 
ternal evolution of the satellites apart from the bound mass 
fraction. 



3.2.1 Satellites with M s = 0.1 M d 

Fig. ^ shows the evolution of the satellite mass for different 
initial orbital inclinations for satellites with M s = M s (0) = 
0.1 Md and eccentricity e ~ 0.7. From this figure we can as- 
sert that: (i) The satellites are disrupted completely at about 
the same time they reach the galactic disc (Fig. |3|) . (ii) For 
all the models, the survival time is, at least, 1 Gyr (25%) 
longer than the equivalent simulations of VW (upper panel 
of Fig. kk). We consider this difference to be indicative of 
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Figure 4 - continued c: As Fig. jlja for satellites with M s 
0.2M^. (Note that the time-axis has changed scale.) 



the uncertainty intrinsic to methods that approximate colli- 
sionless dynamics. The difference comes about, in part, due 
to different numbers of particles, but also due to the spa- 
tial resolution of the method. Prugniel & Combes (1992) 
and Whade & Donner (1996) find that dynamical friction 
is artificially increased due to numerical noise if the particle 
number is small. Similar differences were also noted in the 
computations by Klessen & Kroupa (1998) of satellite ha- 
rassment using different codes. However, we observe that the 
range of disruption times for our models G1S1 (as used by 
VW) is approximately the same, indicating that disc effects 
are well reproduced by our code and giving confidence to 
the following results we obtain using flattened DMHs. (iii) 
Flattened DMHs spread the range of disruption times. In 
Fig. ^a we can see that, for satellites with M s = 0.1 Md em- 
bedded within spherical DMHs this range is ~ 0.9 Gyr (up- 
per panel) , polar satellites having the longest survival time. 
For satellites with the same mass but within flattened DMHs 
the range grows to ~ 1.9 Gyr (lower panel), (iv) Satellites 
with a high orbital inclination within flattened DMHs have 
longer survival times than satellites within spherical DMHs 
with the same initial orbit. For instance, taking the polar 
satellite as the extreme case, G2S190 survives ~ 0.6— 1 Gyr 
longer than G1S190. (v) Satellites with low orbital inclina- 
tion suffer the contrary effect: those within spherical DMHs 
survive longer than those within flattened DMHs. Taking 
the prograde and coplanar orbit as the extreme case, G1S100 
survives ~ 0.4 Gyr longer than G2S100. 

In Fig. |b we compare polar and coplanar satellites 
within flattened and spherical DMHs with orbital eccentric- 
ity e w 0.45 and to obtain an indication of the depen- 
dency of the life-time on e (orbits with intermediate incli- 
nation also have intermediate survival times, Fig. ^a). As 
expected, less eccentric orbits lead to longer survival times, 
since the perigalactic distance is larger and, moreover, tidal 
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Figure 5. The time r when the satellite mass reaches 10 per cent 
of its initial value, M s (t) = 0.1 M s , is plotted vs the initial or- 
bital inclination. Upper panel is for satellite models SI in primary 
galaxies Gl and G2, whereas the lower panel shows the results 
for satellites S2. Note that in all cases r increases with increas- 
ing 8 < 90° for galaxies embedded in a spherical and a flattened 
DMH, due to dynamical friction on the disc. The effect of this is 
particularly nicely seen from the different slopes, dr/d9, for pro- 
grade (6 = - 90°) and retrograde (9 = 90 - 180°) orbits. The 
increase is significantly larger for satellites orbiting in flattened 
DMHs, and becomes larger for decreasing orbital eccentricity (Ta- 
ble |i|) and decreasing satellite mass, which allows longer coupling 
of the satellite to the anisotropic velocity field in the DMH. 



forces are weaker. Furthermore, the survival times show a 
larger spread. Less eccentric orbits survive longer, so that 
anisotropic dynamical friction has a longer time to act. We 
can see that coplanar satellites within a spherical DMH 
(model GISlOOe) survive ~ 0.6 Gyr longer than a coplanar 
satellite within a flattened DMH (model G2S100e), while 
the survival time of a polar satellite within a spherical DMH 
(model GlS190e) is ss 1 Gyr shorter than the correspond- 
ing satellite in the flattened DMH (model G2S190e). Thus, 
the range of survival times increases from about 1.5 Gyr to 
3 Gyr. This range becomes even larger for circular orbits. 

This state of affairs is summarized in Fig. |H] for all 
satellite models, whereas Table ^ compares the decay times 
for SI satellites in dependence of the orbital eccentricity and 
inclination. The table nicely shows that the survival time 
increases significantly with decreasing eccentricity. It also 
shows that oblate DMHs lead to consistently larger differ- 
ences, At, between the decay times for polar and coplanar 
orbits, At consistently being approximately 100 per cent 
larger in flattened DMHs than in spherical DMHs (Ar D bi ~ 
2 AT S ph). This is the key result of this study. 

3.2.2 Satellites with M s = 0.2 M d 

The temporal evolution of satellite masses with M s =0.2 Md 
is shown in Fig. ^c. There are no significant differences in 
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Table 4. Summary of decay times for satellite models SI (M s = 
0.1M|j) in oblate (obi) and spherical (sph) DMHs with different 
initial orbital eccentricity e and orbital inclination 9. tq is the 
decay time when the satellite that is initially on an orbit with 
inclination 9 = 0° has lost 90 per cent of its mass, whereas rgo is 
the decay time for polar orbits (9 = 90°). 



survival times for satellites in spherical and flattened DMHs 
if r a = 55 kpc. At the same time, the dependency on the in- 
clination decreases, causing the range to be narrower in both 
cases. The cause is the fast decay of the satellites, so that the 
anisotropy of the DMH's velocity dispersion does not have 
enough time to act. To better assess this, we introduce a 
set of computations selecting larger initial apo-galactic dis- 
tances (models G3 and G4). The cut-off radius of the Galaxy 
is increased, which changes the rotational curve (see Fig. |l|). 
The results are also plotted in Fig. |^c. A similar spread of 
survival times as for models with M 3 = 0.1 Ma and 'G2' 
flattened DMHs becomes evident; the range of disruption 
times for spherical (G3) and flattened DMHs (G4) are, re- 
spectively, ~ 1 and ~ 2 Gyr. 

The results concerning the disruption times seen on Fig. 
4c between small and large DMHs (G1/G3 and G2/G4 pairs 
displayed on Fig. 4c, bottom panel) are related to one an- 
other as follows. DMHs G3 and G4 have the same mass 
as Gl and G2, but are more extended by a factor r\ = 
133kpc/84kpc = 1.58 (Table [j]). This implies that the dy- 
namical time-scale (oc 1/y/Gp), i.e. the periods of satellites 
on equival ent orbits, ar e longer in haloes G3 and G4 by 
a factor y/ p{G2) / p{Gl) = 1.58 3/2 = 2. Orbits in G3 and 
G4 equivalent to those in Gl and G2, respectively, are or- 
bits with semi-major axes extended by rj in a homologous 
mapping of the systems. Our satellite orbits, however, have 
apo-galactic distances in G3 and G4 twice as large as in 
DMHs Gl and G2. The orbital times of models G3S2nn 
and G4S2nn are in total 1.58 3/2 x 2/1.58 = 2 x 2/1.58 w 2.5 
times longer than models of satellites in DMHs Gl and G2. 
This is approximately what we observe from comparing the 
curves on Fig. 4c with DMHs G1/G3 or G2/G4. 

On the top panel of Fig. ^jc, the time when 
M(*)/M(0) » 0.10 is t w 2 Gyrs for all Gl models. If the ho- 
mologous transformation applied strictly, the curves for the 
G2 halo models should approach 5 Gyrs when M(t)/M(0) = 
1/10. The fact that they are spread between 4 and 5 Gyrs, 
and thus deviate from the homologous map, indicates that 
the disc and bulge, which were left unchanged, play an im- 
portant role in the mass decay rate of the satellites. Further- 
more, the spread in destruction times between models is a 



factor ~ 2, from 1 Gyr (Gl models) to 2 Gyr (G3 models), 
suggesting that the time-scales for orbital decay are con- 
trolled by the DMH, while the combined tidal field of the 
disc and bulge contributes mainly to mass stripping. Sim- 
ilar conclusions would apply for the G2/G4 models shown 
on the bottom panel of the figure. 

3.2.3 Prograde versus retrograde orbits 

Results for models with spherical DMHs may be divided 
into two according to whether the orbit of the satellite is 
aligned with the disc's angular momentum vector (prograde) 
or anti-aligned (retrograde). Keeping the initial satellite ve- 
locity vector unchanged, a prograde orbit is found for an 
initial orbital inclination angle 0° < 8 < 90° , and retrograde 
orbits in the cone 90° < 9 < 180°. 

Table 3 lists four models with spherical Gl DMHs 
and ellipticity e = 0.7 (top segment in the Table). Models 
G1S100 and G1S1180 are respectively prograde and retro- 
grade with respect to the disc, but are otherwise identical. 
From Fig. 4a (top panel) we find for these two simulations 
a 90% mass-loss after w 3 Gyr and 3.5 Gyr, respectively, 
an increase of nearly 20% ; a similar conclusion applies for 
models G1S145 and G1S1135. These findings are qualita- 
tively similar with those of VW: (i) Satellites on prograde 
orbits lose angular momentum faster than their retrograde 
counterparts, leading to more rapid decay, (ii) Polar orbits 
have a similar decay rate as retrograde orbits, as found from 
comparing model G1S190 and G1S1135, Fig 4a. This implies 
that our treatment of the live disc captures the essential 
physics relevant for this work. 

Figure ^ summarizes the findings for decay rates for the 
simulations performed. Point (i) above also applies to flat- 
tened DMHs. However, Fig. ^| suggests in this case that the 
difference in decay rates between prograde and retrograde 
orbits is reduced by about 80 % for flattened DMHs. 

For spherical DMHs the above results can be under- 
stood partially by considering Chandrasekhar's expression 
(Chandrasekhar 1960) for dynamical friction, 

4nG 2 M!(t)p« Vs)]nA 
Fdf = Vs ' ( 5 ) 

Av = \v" s — Vm\ being the relative velocity between the satel- 
lite and the disc particle background, v m is the disc particle 
velocity and p(< v s ) the density calculated only for those 
particles with velocity less than the satellite's velocity v s , 
and InA the Coulomb logarithm, defined as A = p max /p90. 
In this expression, p max is the maximum impact parameter 
(conventionally taken as the half-mass radius of the system), 
and p$o the minimum impact parameter. Since these quan- 
tities are not well defined, the Coulomb logarithm remains, 
to a certain degree, an adjustable parameter. Recent self- 
consistent computations with different iV-body codes leads 
to InA = 1.5 -> 2 (VW; Fellhauer et al. 2000). 

The different decay rate between prograde orbits and 
their retrograde counterparts is caused, in part, by the disc's 
dynamical friction when the satellite is near perigalacticon. 
Retrograde orbits have a much higher relative velocity Av 
due to the disc's rotation and, therefore, they suffer a smaller 
drag force. The bulge or the DMH's dynamical friction make 
no differences since both are non-rotational and spherical, 
which also explains the small differences of decay rates be- 
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tween the polar and the retrograde case (in both cases dy- 
namical friction through the disc can be neglected compared 
to the DMH's dynamical friction). In addition to dynamical 
friction, resonances between the satellite and the disc influ- 
ence the orbital decay, but a detailed analysis goes beyond 
the aim of this work. As for the different decay rates de- 
pending on the satellite's mass, the dynamical friction force 
varies with , so that satellites with M B = 0.2 Md suffer a 
four times larger friction force than those with M s = 0.1 Md- 

3.3 The orbital inclination 9 

Binney (1977) extended the dynamical friction force (|B]) 
to non-isotropic velocity fields. He showed how anisotropic 
friction leads to orbit alignment with the velocity ellipsoid 
plane of symmetry of the host galaxy. Here disc and DMH 
spheroids define a unique z — plane of symmetry, common 
to both mass distribution and velocity ellipsoid. We may, 
therefore, anticipate enhanced satellite orbit alignment rela- 
tively to Binney's analysis, due to the non-uniform, aspher- 
ical mass profile. 

In Fig. |^ we graph the time-evolution of the direction 
angle 9 for a set of simulations with oblate G2 DMHs (q^ — 
0.6) and SI satellites (solid lines on the figure) as well as 
two reference runs with spherical Gl DMHs (dotted lines on 
the figure). 

The average of the orbital inclination 9(t) decreases 
monotonically in time for satellites orbiting in flattened 
DMHs which have initially 9 7^ 0° or 90°. The decrease 
in 9(t) is more appreciable for smaller values of 9(0). This is 
seen for instance by comparing the curves with 9(0) — 15° 
and 30° to the solutions with 9(0) = 60° and 90°. For the 
latter, polar orbit, no decay of 9(t) is observed for the dura- 
tion of the integration, whereas for the 9(0) = 15° case the 
orbit aligns fully with the plane of symmetry of the system 
(coincident with the disc of the host galaxy). 

By contrast, satellites orbiting in spherical DMHs show 
little or no decay of 9(t), for all initial values of 9 (dot- 
ted lines, Fig. ^). This clearly indicates that the anisotropic 
DMH, and not the disc, drives most of the orbital evolution 
and alignment, since in all cases a galactic disc is present. 

The figure also reveals periodic oscillations of 9(t) for 
satellites on inclined orbits, of frequency approximately in 
tune with the satellites' orbital motion. Inspection of the fig- 
ure shows this to be the case for systems with either spherical 
or flattened DMHs. Note that no such oscillations in 9(t) is 
observed for polar or co-planar orbits. We examine the origin 
these oscillations, distinguishing two phases of time around 
t = 2 Gyrs. 

For < t < 2Gyr the satellites orbital radius r s 3> Rd- 
Over this interval of time, the orbits are such that those 
obtained for flattened DMHs lead to much larger oscillations 
in 9(t) compared with the solutions with spherical DMHs. 
We therefore attribute these oscillations to torques from the 
DMH acting on the satellites 

r = rxV * = *fc-^) ^ 

which by symmetry arguments must lie in the plane of the 
axi-symmetric galaxy. The torque V is positive or negative 
according to the phase of the orbit. 

For t > 2Gyr the situation is similar for all calculations, 



independently of the morphology of the DMH. Thus the os- 
cillations we observe clearly for flattened-DMH orbits are 
now noticeable for the solutions with spherical DMHs, too. 
In this phase of evolution, r s ~ Rd or less so that the disc po- 
tential contributes most of the force felt by the satellite and 
hence the torque F acting on it. At this stage, a coupling be- 
tween the disc response and the satellite motion is expected: 
we observed that these oscillations are highly softened in cal- 
culations with a static disc and bulge. Since the orbital angu- 
lar momentum L w r a v s m 3 and AL = Fdt « r 3 GY>(r s /v s ), 
where E is the disc's surface density, both L and the angu- 
lar momentum accrued AL over one revolution will be of 
comparable magnitude if Vg ~ GMd/r s , i.e. when the disc 
potential is the predominant contributor to the force act- 
ing on the satellite. The direction angle 9(t) varies therefore 
wildly towards the end of the simulations in all cases save 
the coplanar 9(0) = 0° one, for which r = at all times. 

The oscillations or periodic fluctuations we have dis- 
cussed are subject to enhancements owing to our choice 
of a grid numerical method of integration. The Cartesian 
grid code limits the vertical resolution of a thin disc. Con- 
sequently the response of the disc to heating by the satellite 
is not correctly quantified. Furthermore, once the remnant 
satellite has merged with the disk, the position of it's cen- 
tre of density becomes ill-defined by virtue of the satellite 
ceasing to exist as a bound entity; 9(t) will reflect this un- 
certainty for 1 5i 2 — 3 Gyr. With 32 mesh points spread over 
a length of 3Rd, the position of the centre-of-density and 
the disc structure are resolved to I ~ 3Rd/32 ~ r 3 /10 when 
r s ss Rd- Hence the error on the angle 9 may be estimated 
to be sin 9 w 9 = l/r s ~ 1/10 or 5° approximately. This 
puts into perspective the magnitude of the oscillations seen 
on Fig. 6 for t <; 2 — 3 Gyr, though without accounting for 
them fully. This leads us to conclude that the physical effect 
of the torque F by the disc on the satellite is qualitatively 
correct, although the quantities somewhat uncertain. 

3.4 Orbital eccentricity 

In Fig. |?] we plot the eccentricity evolution for satellites with 
mass M a = 0.1 Md- The eccentricity is calculated by fitting 
straight lines to r a (t) and r p (t) (e.g. Fig. ^a) and interpo- 
lating e(t) until the satellite has 10 % of its initial mass. 

The orbital eccentricity does not remain constant as 
dynamical friction shrinks the orbit. The evolution of e(t) 
depends on e(t = 0) and 9(t = 0), but from Fig. ^b we 
observe that the general behaviour is for the orbits to be- 
come more radial with time. The only clearly evident excep- 
tion is prograde model G2S100 (e(0) = 0.7), which shows 
a pronounced decrease of e(t). In this case, dynamical fric- 
tion from the flattened DMH plus disc is so large that the 
apo-galactic distance decreases much faster than the peri- 
galactic distance, as we can observe in Fig. [?]a. Close inspec- 
tion shows that this is merely the extreme of a general trend. 
Comparing the co-planar prograde orbits (9 — 0°: GnSlOO, 
GnSlOOe, GnSlOOc; n= 1, 2) with the polar orbits (9 = 90°: 
GnS190, GnS190e, GnS190c), it is evident that the former 
show a stronger sensitivity on initial eccentricity than the 
latter. The effect is such that circular co-planar prograde or- 
bits gain eccentricity rapidly, whereas eccentric (e(0) « 0.7) 
prograde co-planar orbits circularize. Disc-satellite coupling 
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Figure 6. Evolution of the orbital inclination for models G2S100, 
G2S115, G2S130, G2S190 (full lines, satellites within the flattened 
DMH) and G1S145, G1S190 (dotted lines, satellites within the 
spherical DMH), until they retain 10% of their initial mass. 



via dynamical friction and induction of spiral modes in the 
disc and associated transfer of angular momentum between 
satellite and disc are the likely reason, but we do not dwell 
longer on this, as disc-satellite coupling is not the main topic 
of this work, which in any case does not resolve the disc ver- 
tical structure. We merely state here that the data in Fig. 0b 
suggest that there is a stable eccentricity, e sta b w 0.5 for co- 
planar prograde orbits in our flattened DMH, such that e(t) 
increases when e(0) < e s tat>, whereas e(t) slighty decreases 
when e(0) > e sta b, and e(t) « e sta b for all t until satellite 
disruption, if e(0) ~ e sta b- However, if e(0) is smaller than 
e s tab, e(t) does not remain close to e sta b once it has reached 
this critical eccentricity. For our spherical DMH, e sta b ~ 0.6. 

This behaviour agrees with that found by van den Bosch 
et al (1999). They perform numerical calculations using a 
galaxy models similar to Gl, with satellite masses on the 
order of that of our models SI. They observe that the eccen- 
tricity remains remarkably constant for e(0) > 0.6. Unfor- 
tunately, they do not include calculations with lower initial 
eccentricity that we can compare with. We note in passing 
that Prugniel & Combes (1992) already found that initially 
circular orbits rapidly acquire eccentricity, as is also evident 
in the calculations of Fellhauer et al. (2000) . 
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Figure 7. a: Evolution of the apo-galacticon and perigalacticon 
distance (in kpc) for the models with e(t = 0) = 0.7 until the 
satellite has 10 % of its initial mass. 
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Figure 7 - continued b: The eccentricity evolution for the models 
shown in Fig.|^.. e(t) = 1 — fp(t)/r a (t) and is measured by fitting 
the curves r a (t) and r p (t) by straight lines. We also plot it for the 
models with e(0) = 0.45 and e(0) = (circular orbits), calculated 
following the same scheme. 



3.5 Orbital precession 

The orbital plane of a satellite and its unbound particles 
precesses in a flattened potential which smears out the tidal 
debris stream. The precession angle, P(t), is calculated by 
projecting the orbital angular momentum vector onto the 
galactic xy plane and measuring its change with time. In 
Fig. B we plot P for some of our models. The precession, 



dP/dt, increases at later times due to the anisotropy of the 
disc's potential, the satellite having decayed to its vicinity. 

As expected, flattened DMHs lead to larger preces- 
sion. Comparing models CIS 145 (satellite within a spher- 
ical DMH) and G2S145 (satellite within a flattened DMH), 
we observe that the change of P is, respectively, ~ 50° and 
~ 150°, i.e, approximately three times larger at t = 3 Gyr. 
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Figure 8. The precession angle P for some of our models. 

Since the DMH is spherical for models with Gl the pre- 
cession of the orbital plane is due to the disc gravitational 
quadrupole moment. The orbital plane precesses faster the 
smaller its inclination is, orbits with 6 5s 45° precessing by 
180° in 3 Gyr. Polar orbits do not precess at all. 



3.6 Tidal streams 

The accretion history of the Milky Way and other major 
galaxies leaves signatures in the form of old tidal streams in 
the DMHs of these galaxies as found in observational sur- 
veys such as that of Dohm-Palmer et al (2001), or Martinez- 
Delgado et al (2001). The detection of the Sagittarius dwarf 
tails (Iabata et al 1994) therefore likely is a generic features 
of large galaxies. 

Theoretical models of this process have shown good 
agreement with observations (Helmi & White 1999; Zhao 
et al. 1999; Helmi & de Zeeuw 2000). The changes in orbital 
inclination 6 and the orbital precession in flattened systems 
imply that the tidal debris emanating from a disrupting 
satellite will significantly spread out in 8, which will make 
reconstruction of the accretion history of a major galaxy 
difficult if its DMH is flattened. 

In Fig. ^ we plot the deviation angle of the satellite's 
particles from the initial orbital plane in three time snaps. 
This is done for models G1S145 and G2S145 (Fig. 0(0) = 
45°), and for G1S190 and G2S190 (Fig. §d, 0(0) = 90°). The 
first time-snap shows satellite particles after first passage 
through perigalacticon at t = 0.62 Gyr, the second one is 
at an intermediate time (t = 1.52 Gyr) while the last frame 
is at a late stage of the satellite orbit. The debris does not 
remain in the initial orbital plane. This effect becomes more 
pronounced the closer the satellite is to the galaxy's cen- 
tre, when the mass loss (Fig. ^) and the oscillations of the 
orbital inclination (Fig. p|) primarily occur, and the larger 
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Figure 9. a: Deviation angles for all satellite particles from the 
initial orbital plane (8 = 45°). The left column depicts model 
G1S145 (spherical DMH), and the right column shows G2S145 
(flattened DMH). Rows show three time snaps (given in Gyr). In 
the last one, the satellite has been fully destroyed. 

the number of perigalacticon passages is. From Fig. ^a we 
also observe that the deviations from the orbital plane are 
enhanced when the DMH is flattened since satellite orbits 
within oblate DMHs align with the symmetry plane (i.e. 
6(t) — > 0). Fig. |^b shows that the spread of satellite debris 
is much smaller for satellites in polar orbits than for those 
with intermediate inclinations, since inclination decay and 
oscillations vanish for polar orbits. 



4 CONCLUSIONS 

In order to assess the importance of dynamical friction in 
extended oblate DMHs on the distribution of satellite galax- 
ies around their primary, we perform self-consistent iV-body 
computations of satellite galaxies with masses amounting 
from 10 to 20 per cent of the primary's disc. The satel- 
lites are placed on different orbits in spherical and flattened 
DMHs that have embedded galactic discs and bulges. 

The calculations with spherical DMHs lead to results in 
good agreement with those obtained by Velazquez & White 
(1999). Modest differences in quantities are attibuted to the 
increased mass resolution of our calculations compared with 
theirs, as well as different linear resolution (grid size versus 
smoothing length of their TREE algorithm). 

Satellites evolving in spherical DMHs on prograde orbits 
relatively to the primary galaxy's disc rotation decay faster 
than satellites on retrograde orbits or on polar orbits. This 
results from orbital resonances between the disc and the 
satellites. 

Of particular interest, however, is that our results 
demonstrate that non-isotropic dynamical friction in flat- 
tened DMHs works as a removal mechanism of satellites with 
low-inclination orbits, whereas it enhances the survival time 
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Figure 9 — continued b: As Fig. but for models G1S190 (spher- 
ical DMH) and G2S190 (flattened DMH), with initial inclination 
6 = 90° . 



of satellites on near polar orbits. Thus, satellites on polar 
orbits survive about 70 per cent longer than satellites on 
orbits that have a small inclination relative to the primary 
galaxy's disc (Table irrespective of the relative orbital 
sense (Fig. in an oblate DMH with axis ratio = 0.6. 
This is the key result of this investigation. 

This result helps understand the distribution of dwarfs 
galaxies in the Milky Way. Since they are mainly distributed 
near the galactic pole (Carney et al. 1987) we may infer a 
selection of survivor dwarfs from a primordial population. 
The accelerated orbital decay and alignment with the disc 
of dwarfs within a flattened halo would go some way towards 
accounting for the data. However if the masses deduced for 
these satellites ( 10 s solar, compared with 10 9 for our mod- 
els) is a good measure of their mass at the formation time, 
our computations indicate times as long as a Hubble time 
for effective mergers. Discrepancies in timescale may well be 
accounted for if we substitute for the isothermal halo the 
more concentrated NFW (Navarro, Frenk & White 1995) 
models or haloes with a steeper cusp (Moore et al. 1998): 
when each halo model is scaled to the same integrated mass 
inside the solar radius, the particle velocity dispersion in 
these models drops faster with radius than for isothermal 
spheres. Because of the strong dependence of friction on ve- 
locity dispersion, this would reduce the timescale for orbital 
decay very much and offset the effect of reduced satellite 
masses. We have not, however, performed calculations with 
different halo mass profiles. 

Our computations further show that satellites on or- 
bits with eccentricity e ^ 0.5 and with masses larger than 
10 per cent of their primary galaxy's disc merge within only 
a few Gyr with the primary galaxy. The time it takes to 
merge increases with decreasing orbital eccentricity (Fig. [5J). 
We therefore deduce that massive satellites around distant 
galaxies, such as typically enter the samples that show the 



Holmberg effect, may be preferentially on near-circular po- 
lar orbits or on orbits with apo-galactica further away from 
their primary galaxy than about 130 kpc. 

The calculations also suggest that there exists a critical 
eccentricity, e sta b, for co-planar prograde orbits such that 
initial eccentricities that are close to e a tab remain within 
the vicinity of this 'stable' value, whereas initial orbital ec- 
centricities that differ from e s t a b evolve towards the critical 



value but the orbit keeps evolving past e sta b (Section 3.4). 
For our oblate DMH, e ata b ~ 0.5, whereas for the spherical 
DMH e stab w 0.6. 

We also note that the high precession rates of satellite 
orbits in flattened DMHs and the decrease in orbital inclina- 
tion leads to tidal debris streams being completely smeared 
apart for initially inclined orbits. 

We want to comment that, despite our use of only 
two values for the satellite mass in our calculations, this 
range seems to be representative to reproduce the typical 
mass of the satellite that Holmberg (1969) and Zaritsky 
& Gonzalez (1999) find in their observations when the ini- 
tial apo-galactic distances is selected properly (Ibata et.al 
2001). As Tormen (1997) finds in his numerical calculations 
of hierarchical galaxy clusters history, more massive satel- 
lites (~ 10 11 ) are unlikely to survive due to the large drag 
force they suffer. On the other hand, though less massive 
satellites (~ 10 s solar) feel a negligible drag force, they are 
quickly disrupted after some peri-galacticon passages due to 
their low binding energy. 

This paper has sought to quantify the effect of aspher- 
ical DMHs on the orbits of galactic satellites. The analysis 
suggests enhanced Holmberg decay, yet what can we say of 
a population of satellites as a whole? Our model satellites 
require a few orbits around the host galaxy if dynamical 
friction is to be effective. Thus within one Hubble time a 
satellite would require = 5 revolutions (say) or t = 2 Gyr 
for a single revolution at most. In the Milky Way the orbital 
time t — 200 Myr at r = lOkpc; assuming an isothermal 
halo with p tx r~ 2 , the critical orbital time t = 2 Gyr would 
be found at r = 50 kpc or so. In other words, satellites that 
are too far from the host galaxy will not have time to ex- 
perience dynamical friction and hence will not have suffered 
Holmberg decay. On the other hand, satellites closer to their 
host galaxy will merge quickly through the process described 
here. Zaritsky et al. (1999) have noted that satellite popu- 
lations tend to remain isotropically distributed for satellites 
with r > 50 Kpc. 

A more elaborate study is under way, and ultimately 
we aim at making a statistical study of a modelled observa- 
tional sample to infer if the Holmberg effect can indeed be 
produced by flattened DMHs. 
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